ARPA  ORDER  NO.  189-61 


.  'Memorandum 

^RM-4809-ARPA 

^  DECEMBER  1965 


*  i 

■ 


DETECTION  OP  SONAR  SINUSOIDS  OF 
UNKNOWN  FREQUENCY  AND  KNOWN 

OR  UNKNOWN  PHASE 

F.  S.  Hill,  Jr. 


^  C  I  E  A  R  !  r,  H  n  ?|  t:  r 
r’OH  FEDEpau  ^  . . 

Eardccp*  ~  u  •  •  _ - 

75'  ,1^ 


^ '  R  » ■  I  »  •  • 

w  a  -1>- , , 


U  L» 


PREPARED  FOR-  / 

ADVANCED  RESEARCH  PROJECTS  AGENCY 


lie 


RUnD 


SANTA  MONICA  •  CAIIFOUNIA- 


I 


ARPA  ORDER  NO.  189  01 


MEMORANDUM 

RM-4809-ARPA 

DECEMBER  1965 


DETECTION  OP  SONAR  SINUSOIDS  OP 
UNKNOWN  PREQUENCY  AND  KNOWN 

OR  UNKNOWN  PHASE 

F.  S.  Hi?l,  Jr. 


This  research  is  supported  by  the  Advanced  Research  Projects  Agency  under  Contract 
No.  SD-79.  Any  views  or  conclusions  containi.d  in  this  Memorandum  should  not  be 
interpreted  as  representing  the  official  opinion  or  policy  of  ARPA. 


DISTRIBUTION  STATEMENT 
Distribution  of  this  document  is  unlimited. 


7^ 


m]d 


<  700  S  ’ 


)ANTA  MONiCA 


'AllfOtNl*  '  toaD** 


Approvecj  for  release  by  the 
Federal  Scientific  and  Techni 


CiearinBiiouse  for 
cal  !nfoffn;..ion 


iii 


PREFACE 

In  conjunction  with  RAND's  study  of  Defense  Against  Submarine- 
Launched  Ballistic  Missiles  for  the  Advanced  Research  Projects  Agency, 
background  investigations  of  theoretical  methods  for  calculating  the 
performance  of  nondirectional  passive  sonobuoys  are  being  conducted. 
The  particular  investigation  presented  in  this  Memorandum  evaluates 
the  performance  of  a  theoretically  optimum  processor  for  detecting 
a  sine  wave  of  unknown  frequency  and  unknown  phase  in  gaussian  noise. 

This  study  should  be  of  interest  to  analysts  in  the  fields  of 
sonar  and  radar  as  well  as  to  researchers  in  the  general  field  of 
signal  detection  theory. 
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SUGARY 

This  Memorandum  considers  the  problem  of  detecting  a  constant 
sine  wave  of  unknown  frequency  and  amplitude  in  gaussian  noise.  In 
particular,  it  is  assumed  that  the  sinusoid  may  appear  at  any  one  of 
a  finite  number  of  known  frequencies,  and  the  probability  of  its 
'ccurrence  at  each  of  these  frequencies  is  assumed  to  be  equal.  Two 
cases  are  treated  here.  The  first  assumes  that  although  the  frequency 
is  not  known,  the  phase  of  the  signal  is  known,  thus  allowing  co¬ 
herent  detection.  The  second  acknowledges  that  the  initial  phase 
indeed  could  not  be  known,  and  an  analysis  of  the  incoherent  detector 
is  made.  It  is  shown  that  for  large  output  signal- to -noise  ratios, 
the  problem  in  both  cases  becomes  that  of  detecting  one  of  m  approxi¬ 
mately  orthogonal  signals  in  a  noise  background*  The  magnitude  of 
the  error  in  the  orthogonality  approximation  is  considered.  A 
physical  realization  of  an  approximately  optimum  detector  structure 
is  studied  in  some  detail,  and  the  effect  of  finite  observation  time 
is  considered. 

The  results  indicate  that  the  difference  between  the  two  cases 
is  small  and  quite  predictable.  Thus,  in  studies  of  this  kind  it  can 
b(!  assumed  that  the  initial  phase  is  known,  the  gaussian  character  of 
the  quantities  that  arise  may  be  retained,  and  finally,  the  answers 
can  be  adjusted  to  account  for  the  actual  lack  of  knowledge  concerning 


the  phase. 
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I,  INTRODUCTION 

The  problem  considered  here  is  that  of  detecting  the  presence 
of  a  submarine  target  by  a  single,  unattended,  omnidirectional  hydro¬ 
phone.  It  is  presumed  that  such  a  target  emits  both  wideband  noise 
and  a  set  of  sine  waves  of  different  frequency.  This  Memorandum 
considers  the  detection  of  the  target  by  means  of  the  sine  waves. 

Since  the  exact  frequency  of  each  sinusoid  will  not  be  known  in 
advance,  tue  problem  becomes  that  of  detecting  one  or  '^ore  sinusoids 
of  unknown  frequency  in  a  noise  background. 

As  a  first  step  toward  a  solution  of  this  problem,  the  question 
of  detecting  a  single  sinusoid  of  unknown  amplitude  and  frequency  in 
a  background  of  gaucsian  noise  is  considered.  The  noise  is  assumed 
to  have  a  known  spectral  density  and  is  assumed  stationary,  at  least 
for  the  duration  of  each  decision  period.  In  practice,  the  frequency 
of  the  sinusoid  may  only  be  known  within  a  band  of  frequencies,  and 
may,  with  equal  probability,  appear  anywhere  within  this  band.  As 
a  useful  approximation  tc  this  situation,  the  assumption  is  made  here 
that  the  frequency  of  the  sinusoid  has  a  discrete  distribution;  that 
is,  it  can  appear  at  one  of  a  finite  number,  m,  of  possible  frequencies. 
If  this  assumed  number  m.  is  large,  then  the  assumption  is  not  greatly 
different  from  that  of  assimiing  a  continuous  distribution  of  fre¬ 
quency  over  the  band  of  uncertainty. 

It  is  seen  later  in  this  Memorandum  that  the  false-alarm  rate, 
b'iing  determined  only  by  the  noise  statistics,  is  independent  of 
signal  amplitude,  and  that  detectability  of  the  target  increases 
monotonically  with  signal  amplitude.  Thus,  if  signal  amplitude  is 
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not  known,  but  is  known  to  remain  steady  over  the  length  of  a  de¬ 
tection  period,  then  the  system  threshold  is  set  by  the  allowable 
false-alarm  rate,  and  its  detection  performance  will  depend  only  on 
the  actual  signal  amplitude.  Hence,  a  constant  signal  amplitude  of 
unknown  level  can  be  assumed,  as  done  in  this  Memorandum,  giving 
final  answers  in  terms  of  the  actual  signal  amplitude.  The  signal 
amplitude  may,  of  course,  -e  converted  into  target  range  if  the 
signal  level  emanating  from  the  target,  the  transmission  loss  versus 
distance  and  the  noise  level  at  the  receiver  are  known. 

The  first  part  of  the  Memorandum  deals  with  the  detection  of  a 

sinusoid  of  unknown  frequency  but  known  phase  in  gaussian  noise.  The 

known  phase  assumption  is  clearly  unrealistic  if  the  frequency  is 

not  known,  but  the  mathematics  are  tractable,  and  it  is  known  from 

studies  of  knovm- frequency  sine  wave  detection  that  the  difference 

* 

between  coherent  and  incoherent  techniques  is  small  and  predictable. 
The  second  part  demonstrates  this,  dealing  with  the  detection  of 
sinusoids  of  both  unknown  frequency  and  phase.  This  study  concludes 
that  if  work  is  pursued  knowing  the  phase  of  the  sinusoid,  the  effect 
of  the  actual  ignorance  of  phase  may  then  be  included  by  means  uf 
some  loss  factor  which  is  constant  over  a  large  range  of  situationc. 


Confusion  sometimes  arises  concerning  the  terms  "coherent"  and 
"incoherent,"  because  their  roles  in  disrussions  of  active  radar 
systems  are  quite  different  from  their  usage  here.  In  radar  nomen¬ 
clature,  "coherent"  and  "incoherent"  refer  to  methods  of  integrating 
sequences  of  pulses,  either  retaining  or  destroying  the  phase  struc¬ 
ture  between  different  pulses.  In  the  present  context,  however,  there 
are  no  pulses  at  all.  Here,  "coherent"  specifies  actual  knowledge  of 
signal  phase  relative  to  some  reference,  for  all  time.  "Incoherent" 
implies  the  absence  of  this  knowledge,  which  sugge.>ts  the  use  of  a 
statistical  description. 


II.  STATEMENT  OF  THE  PROBLEM 


A  single  hydrophone  monitors  underwater  acoustical  energy  for 
an  observation  time  of  T  sec.  If  a  signal  is  present,  it  is  assumed 
to  have  the  following  form 

s(t)  =  P  cos(ui>t  +  cp)  (1) 

where  P  is  assumed  constant,  the  phase  may  or  may  lOt  be  known,  and 
the  frequency  =  2r\  f^  is  a  random  variable  having  a  discrete 
probability  density  function 

m 

p(a')  =  ^  P(uj^)  6(u)  -  UD^)  (2) 

i=l 

r  1  i  f  X  “  0 

where  6(x)  =  i  n  -r  an  is  the  Kronecker  delta  function. 

'  lOifxfO 

That  is,  the  signal  appears  at  one  of  the  frequencies  and  is 
assumed  to  remain  at  that  frequency  for  the  duration  of  the  obser¬ 
vation  interval.  The  frequency  nay  possibly  "wobble"  slightly  from 
interval  to  interval  but  is  assumed  to  change  so  slowly  that  it  is 
essentially  constant  for  the  period  T. 

The  m  possible  frequencies  at  which  the  signal  may  appear  are 
assumed  to  be  contained  in  a  low-frequency  band  of  total  width  Q. 

The  noise  process  received  by  the  hydrophone,  in  addition  to  the 
signal  (if  present),  is  assumed  to  be  a  stationary,  gaussian,  random 
process  with  zero  mean,  and  has  a  known  spectral  density  N(f) .  The 
detector  will,  be  concerned  with  the  values  of  f  lying  in  (-H,  Q)  ,  and 
thus  N(f)  is  considered  as  being  non-zero  only  in  this  band.  The 
bandwidth  is  then  Q  cps. 
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It  is  known  from  sampling  theory  that  a  waveform  of  duration  T  sec 
and  bandwidth  Q  cps  may  be  approximately  represented  by  a  set  of 
20  T  numbers.  If  the  frequency  band  of  the  wave  extends  to  zero  cps, 
then  these  numbers  are  simply  equally  spaced  amplitude  samples  of  the 
wave.  In  the  following,  it  will  be  mathematically  convenient  to 
treat  the  received  waveform  in  terms  of  these  20  T  amplitude  samples. 
The  actual  detection  system  will  not  perform  any  sampling,  however; 
the  discrete  representation  is  merely  a  technique  to  facilitate  manipu¬ 
lation.  The  information  thus  presented  to  the  detector  is  a  vector  v 
of  n  voltage  samples  where  n  =  20  T.  The  column  vector  v  has  the 
transpose 


=  [v(tp,  vCtp,  ...»  vCt^)]  (3) 

v(t^)  =  n(tj^)  if  signal  is  absent 

v<'t^)  =  n(t^)  +  P  cos(a)t^  +  if  signal  is 

present  at  frequency  u) 

It  is  well  known  that  the  optimum  detector  for  a  signal  con- 
taminaced  by  a  random  noise  is  a  likelihood  ratio  detector. If 
any  of  the  input  quantities  contain  unknown  parameters,  the  likeli¬ 
hood  ratio  must  be  averaged  over  all  possible  values  of  these  quanti¬ 
ties  according  to  their  probability  distribution.  For  instance,  if 

the  signal  contains  the  paramerers  9  and  0  which  are  statistically 

I  z 

independent  but  unknown  a  priori  to  the  observer,  the  detector  must 
form  the  averaged  likelihood  ratio 

=  J  JdPcep  l(v.  e^,  ep 


(4) 
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where  the  quantity  P(6j^)  is  the  cumulative  probability  distribution 
function  of  9^,  and  similarly  for  P(02)*  l(v,  9^^,  82'  likeli¬ 

hood  ratio,  depending  on  the  received  vector  v,  and  assuming  that  the 
signal  parameters  have  the  particular  values  9^^  and  02*  The  proba¬ 
bility  structure  of  the  unknown  parameters  may  or  may  not  be  known 
in  a  particular  case.  In  the  case  where  9^^  is  the  phase  angle  of  a 
sinusoidal  signal,  it  is  reasonable  to  assume  that  the  angle  cp  is 
uniformly  distributed  in  the  interval  0,  2tt.  This  assumption  will 
be  made  in  the  following  when  the  phase  is  assumed  unknown. 

Hence,  in  the  two  cases  to  be  considered  in  this  Memorandum;  the 
coherent  case  where  the  phase  angle  is  assumed  known,  and  the  in¬ 
coherent  case  where  it  is  assumed  to  be  a  uniformly  distributed  random 
variable,  the  two  quantities  must  be  formed 

m 

Coherent  j^(v)  =  S  L(v,  u).)  p(u).)  (5) 

i=l  ^  ^ 

m 

Incoherent  £(v)  =  ^  p((D^)  J  fn  '“i’ 

i=l  ^ 


The  sunmation  performs  integration  over  the  unknown  but  discretely 
distributed  frequency  parameter,  according  to  Eq.  (2).  The  two  cases 
will  be  considered  separately. 

Although  signal  level  is  not  known,  it  is  not  considered  as  a 
random  variable  here,  and  no  averaging  is  performed  over  this  parameter 
because  the  operations  which  the  optimum  detector  performs  are  not 
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influenced  by  signal  amplitude.  The  resulting  detection  performance 
depends  very  strongly  on  signal  strength,  but  the  detector  can  be 
built  with  no  knowledge  of  how  strong  the  signals  will  be.^^^  Noise 
level,  on  the  other  hand,  must  be  known  in  order  .  '  set  a  false- 
alarm  rate. 
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III.  DETECTION  OF  A  SINE  WAVE  OF  UNKNOWN  FREQUENCY  AND  KNOWN  PHASE 

The  case  of  detecting  a  single  sine  wave  of  known  initial  phase 
but  unknown  frequency  has  been  considered  in  some  detail  by  Levesque, 
and  this  section  reports  his  findings.  It  is  realized  that  the 
assumption  of  known  phase  is  unrealistic  when  the  frequency  is  un- 
knoT-m,  but  it  is  believed  that  the  following  analysis  is  useful  for 
two  reasons:  (1)  it  manipulates  gaussian  statistics  and  is  therefore 
tractable,  and  (2)  it  places  an  upper  bound  on  detectability  which 
is  later  (Section  IV)  seen  to  be  easily  related  to  the  unknown  phase 
case  results. 

The  optimum  detector  forms  the  quantity  given  in  Eq.  (5)  and 
compares  it  with  a  preset  threshold.  If  exceeds  this  threshold, 

the  decision  is  made  that  a  signal  is  present.  If  does  not 

exceed  the  threr.hold,  the  signal  is  dismissed.  Levesque  introduces 
the  concept  of  the  decision  hyperplane,  which  permits  a  geometrical 
interpretation  of  the  decision  process.  Considering  the  simple  case 
where  the  signal  may  appear  at  one  of  only  two  possible  frequencies, 
m  =  2,  the  averaged  likelihood  ratio  takes  the  form 

£(v)  =  p(u)^)L^  +  P^^2^^2 
or 

=  p(aj^)  exp(log  L^)  +  pCa;^)  exp(log  L^)  (8) 

where  L^^  =  l(v,  o)^)  and  L^  =  l(v,  for  convenience.  The  form  of 
Eq.  (8)  is  used  as  it  is  usually  more  convenient  to  determine  the 
quantity  log  L^  than  the  likelihood  ratio  itself.  The  optimum 


(2) 
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detector  compares  the  quantity  of  Eq.  (8)  with  a  threshold  k,  and  a 
target-present  decision  is  made  depending  on  whether  -C(v)  exceeds 
this  threshold  or  not.  Thus,  the  "decision  line"  that  separates  one 
decision  from  the  other  is  that  which  satisfies  the  equation 

/(v)  =  k  (9) 

Assuming  that  the  two  possible  frequencies  are  equally  likely, 
p(u)j^)  “  ~  ^9®*  (8)  and  (9)  yield 

log  =  log[2k  -  exp(log  Lj^)  ]  (10) 

as  the  equation  of  the  "decision  line."  This  line  is  shown  in  Fig.  1 
for  various  values  of  the  threshold  k.  As  one  of  the  quantities 
log  or  log  be  omes  small,  the  decision  curve  rapidly  approaches 
an  asymptote  parallel  to  one  of  the  coordinate  axes.  In  an  actual 
system,  a  preset  value  of  k  is  used,  and  thus  one  of  the  curves  of 
Fig.  1  will  form  the  boundary  between  the  signal-accept  and  signal- 
dismiss  regions.  That  is,  the  processed  vector  v  yields  the  quantities 
log  L, ,  log  which  correspond  to  a  point  on  the  graph  of  Fig.  1. 

If  this  pcint  lies  above  or  to  the  right  of  the  decision  line  given 
by  Eq.  (10),  a  signal-accept  decision  is  made.  If  log  or  log  L2 
lies  below  or  to  the  left  of  this  line,  the  signal-dismiss  decision 
is  made.  If  the  number  m  of  possible  frequencies  were  greater  than 
two,  then  the  likelihood  ratio  of  Eq.  (5)  would  be  compared  to  some 
other  threshold.  The  resulting  decision  surface  of  m  dimensions 
would  divide  an  m-dimensional  hyperspace  into  two  regions:  signal- 
accept  and  signal-dismiss.  The  processed  data  would  yield  m  numbers, 
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Fig.  1— Coherent  sinusoid:  decision  curve  for  case  of  two  possible 

frequencies  (m  =  2) 
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log  Lt ,  log  L_ ,  . . . ,  log  L  ,  corresponding  to  some  point  in  the  m- 
dimensional  hyperspace;  and  the  signal-accept  decision  is  made  if 
the  point  lies  in  the  appropriate  region.  Otherwise,  the  signal- 
dismiss  decision  is  made. 

FORM  OF  THE  LIKELIHOOD  RATIO 

The  likelihood  ratio  L(v,  may  be  written  explicitly.  It  is 
well  known  that  for  a  signal  s^(t)  =  P  cos(u.)^t  +  cp)  imbedded  in 
gaussian  noise  with  covariance  matrix  K,  the  likelihood  ratio  has  the 
form 

L(v,  (u^)  =  exp  1^-  K  ]  (11) 

Here  ^  is  the  row  vector  of  time  samples  of  the  sinusoid  at  angular 
frequency  u)^.  The  noise  covariance  matrix  K  is  related  to  the  noise 
spectral  density  N(f)  by  the  relation  [Ref.  3,  p.  103] 

00 

R(t)  =  J  N(f)  exp(juJT)  df  (12) 


K  = 


i.  \  1 

where  ^  ~  between  samples. 
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The  logarithm  of  the  likelihood  ratio  of  Eq.  (11)  has  the  form 

log  L.  =  -  s'  k"^  s.  +  s'  k"^  V  (13) 

The  second  term  indicates  that  the  detector  must  correlate  the  re¬ 
ceived  waveform  with  the  signal  at  the  frequency.  For  example, 

if  the  noise  is  assumed  to  have  a  flat  spectral  density  of  level  N  /2 

o 

2 

volts  /cps  over  the  band  of  all  possible  frequencies  (“Q,  Q) ,  then 
the  matrix  K  is  given  by  N^I  0,  where  I  is  the  unit  matrix,  the  noise 
variance  N  =  and  the  correlation  portion  of  Eq.  (13)  becomes 

n 

s'  K  '  V  =  ^  P  cos(oj^t  +  ip)  v(t  )  (14) 

2P 

=  JT  J  cos(uj^t  +  fp)  dt  (15) 

o  “o 

Equation  (15)  approximates  the  sum  of  Eq.  (14)  by  an  integral  and 
makes  use  of  the  relation  dt  =  At  =  1/(2Q). 

The  first  part  of  Eq.  (13)  may  likewise  be  analyzed,  and  it  is 
seen  to  become 


I? 


n 


-  7  s'  K-l  ^  i 


(16) 


j=l 


1  2  1*^  2  2 

=  -  2  ^  J  P  cos  ((i)^t  +  cp)  dt 


o  o 


2  N 

o 


(17) 


2  jT 

as  long  as  T  is  large  in  comparison  to  each  sine  wave  period, 

It  is  convenient  to  define  the  "detection  index,"  d,  which  satisfies 


d^  -  P^/N^  (18) 

2 

Thus,  d  is  the  signal-to-noise  ratio  at  the  output  of  the  detector 

at  the  end  of  the  observation  period  of  T  sec  and  may  be  used  as  a 

figure  of  merit  in  evaluating  detectability.  .is  is  not  the  signal- 

to-noise  power  ratio  at  the  input  to  the  detector,  which  has  the 
2  . 

value  P  /2N,  but  it  is  the  ratio  after  the  detector  has  processed 
the  received  data.  By  comparing  Input  and  output  signal-to-noise 
ratios,  it  is  seen  that  the  "processing  gain"  is  2Q  T,  the  number  of 
degrees  of  freedom  in  the  sampling  representation.  Note  that  although 
P  and  therefore  d  are  not  known  a  priori,  they  can  be  used  here  as 
if  they  were  known,  and  the  results  can  be  presented  in  terms  of  d 
as  a  parameter. 


PEPJORMANCE  OF  THE  DETECTTOR 


It  remains  to  calculate  the  perfvjrmance  of  the  detector.  This 
requires  findiup  the  probabilities  of  false  alarm  and  detection  for 
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particular  settings  of  the  threshold  and  for  particular  values  of  input 
signal- to-noise  ratio.  In  order  to  do  this,  the  conditional  joint 
distributions  of  the  test  quantities  log  ,  log  must  be  computed 
under  the  three  hypotheses 

no  signal  present,  v(t)  =  n(t)  (19) 

signal  present  at  ud^^,  v(t)  =  n(t)  +  Sj^(t)  (20-a) 

signal  present  at  ~  "*■  (20-b) 

Levesque  has  done  this ,  and  the  derivations  appear  in  Appendix  A. 

It  is  shown  there  that  the  resultant  joint  distribution  is  gaussian 
and  that  the  quantities  log  and  log  are  approximately  indepen¬ 
dent  since  there  is  almost  no  correlation  between  sine  waves  of  dif¬ 
ferent  frequencies  over  a  long  observation  time,  a  time  such  that 
rZC  >  >  1  and  |  T  >  >  1.  This  is  exactly  true  if  ju;^  -  uu^  |  is 

a  multiple  of  2n/T  cps.  The  results  are 


.2  >.2 


2  ^2 


f(x, ,X-/H  )  = 


/  "  9 

^  °  2nd 


2d^ 


(21) 


.2  ^2 


2  ^2 


f(x  ,x  /H.)  =  - exp 

^  2  ^  2Tld‘ 


r  <1  t  ^  ^  <1  ^ 


2d 


(22) 


f(x.  ,x  /H  )  =  - -  exp  I 

^  2  "  2nd^  i 


r  +  d^ 

J  *{^2--  J  I 

.2 


2d 


(23) 


where  for  convenience  x^  ®  log  L^;  i  =*  1,  2. 
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It  is  noted  that  the  shift  in  the  mean  value  of,  say,  , 

between  the  cases  of  signal  absent  and  signal  present,  is  equal  to 
2 

d  .  The  standard  deviation  of  x,  is  d.  Thus,  the  ratio  of  the  shift 

1 

in  the  mean  over  the  standard  deviation  is  given  by  d  itself.  The 

square  of  such  a  ratio  is  frequently  used  as  the  output  signal-to- 

noise  ratio  in  detection  prollems.  This  reinforces  the  choice  of 
2 

d  or  d  as  a  figure  of  merit. 

Using  these  probability  density  functions,  it  is  convenient  to 
use  the  geometrical  approach  of  the  decision  space  to  calculate  the 
performance  probabilities.  These  three  functions  are  represented 
in  Fig.  2,  where  the  contours  of  constant  probability  are  centered 
around  the  mean  values  of  the  density  functions.  For  the  value  of 
threshold  shown,  it  is  seen  that  under  the  hypothesis  noise  only 
is  present,  and  there  is  some  probability  that  the  point  log  L^, 
log  =  Xj^,  X2  will  actually  fall  above  the  threshold  decision 
curve.  This  false-alarm  probability  may  be  found  by  integrating 
the  density  function  of  Eq.  (21)  o',  or  the  entire  region  above  and  to 
the  right  of  the  decision  curve.  Using  Eq.  (10)  of  the  decision 
curve,  this  gives 


=r 

log  2k 


+ 


log  2k  2 

J  dx  — ^  exp  f-  X  +  ] 

ia.  2  /2n  d  -  2d^  ^  ^  ^ 


(24) 


r 


dx. 


exp 


logL2k  -  exp(x  )] 


y2n  d 


— 


1 
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Fig. 2— Location  of  probability  density  functions  relative  to 
decision  regions  under  three  hypotheses  for  m  =  2  case 
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Unfortunately,  the  lower  limit  detemined  by  the  decision  curve  in 
the  above  integ  al  renders  the  integral  unsolvable  in  closed  form. 

A  useful  approximation  to  the  value  of  the  integral  may  be  made  by 
using  the  asymptote  to  the  decisic  curve  instead  of  the  exact  curve 
itself.  To  understand  the  effect  of  this  approximation,  consider 
the  following  argument. 

It  is  seen  from  Fig.  2  that  the  probability  density  functions  of 

2 

Xj^,  x^  are  centered  about  the  values  ±  d  /2.  Thus,  as  the  signal- 

to-noise  ratio  increases,  the  centers  of  these  distributions  move 

away  from  the  origin  and  move  further  away  from  the  decision  line. 

If  ‘'he  distributions  are  far  from  this  line,  its  exact  shape  is  no 

longer  important,  and  the  asymptotes  may  be  used.  Another  way  of 

considering  this  is  to  examine  the  area.  Region  II,  involved  in  the 

approximation.  By  using  the  asymptotes.  Region  II  is  effectively 

added  to  the  "signal  absent"  decision  portion,  when  it  really  belongs 

to  the  "signal  present"  decision  portion.  The  error  then  is  equal  to 

the  probability  density  integrated  over  Region  II.  But,  as  the 

2 

distributions  move  further  away  with  increasing  d  ,  their  value  within 
this  area  decreases  rapidly.  Therefore,  it  may  be  said  in  general 
that  for  large  output  signal-to-noise  ratios,  the  asymptotes  may  give 
a  very  good  approximation.  Since  the  asymptotes  are  parallel  to  the 
coordinate  axes  in  both  this  and  the  m  >  2  cases,  the  use  of  the 
asymptotes  for  the  decision  curve  is  analogous  to  making  a  decision 
based  on  each  quantity  for  log  independently.  Levesque  calls  this 
behavior  the  "band-splitting"  detector,  since  it  treats  the  processed 
received  signal  at  each  of  the  m  different  frequencies  independently. 
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The  detector  approaches  a  maximum- likelihood  detector  for  large  output 
signal-tc-noi 3e  ratios.  Again,  It  should  be  pointed  out  that  the 
signal-to-noise  ratio  in  question  is  that  at  the  output  after  pro¬ 
cessing.  For  small  input  signal-to-noise  ratios,  the  output  signal- 
to-noise  ratio  may  be  increased  many  times  by  a  long  integration  time, 
as  incicated  by  Eq.  (18),  Since  only  in  the  case  of  fairly  large 
output  signal-to-noise  ratios  is  there  any  chance  of  detection,  these 
are  the  primary  cases  of  interest.  Hence,  this  aoproximation  will  be 
used  in  the  following. 

Returning  to  the  evaluation  of  the  false-alarm  probability,  or, 
and  using  the  value  of  the  asymptote  in  the  lower  limit  of  the  inte¬ 
gration  of  Eq.  (24),  log[2k  -  exp(x2)l  =  log  (2k),  a  may  be  evaluated. 
With  a  change  of  variable,  Eq.  (24)  becomes 


where 


.  1 

Ot  -  - 

I*®  2 

[  exp  C 

)  dv 

v/2W‘' 

a 

^  .2 

"v  f*®  ^  2  .. 

f 

••00 

)<lvj  exp(^) 
a 

=  ^  [^1  -  $  (a)  ]  1^3  +  $  (a) 


du 


(25) 


(26) 


a  =  I  +  ^  log  2k 


(27) 


2 


1  ^ 

and  where  $(x)  =  -  e  *'  dt  is  the  normal  probability  integral 

/in  -X 
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(4) 

given  ii  tables.  '  Using  this  resuli;,  the  proper  threshold  nay  be 
set  for  an  allowable  false-alarm  probability  per  observation  time. 

The  probability  of  detection  is  the  probability  that  the  values 
log  L, ,  log  will  lie  above  the  threshold  curve  under  either  of 
the  hypotheses  H2  (for  equally  likely  signals).  This  is  more 
conveniently  found  by  computing  the  miss  probability,  8,  which  is 
equal  to  one  minus  the  detection  probability  3=1-  P^,  the 
probability  that  the  test  quantities  yield  a  point  that  lies  below 
the  threshold,  even  though  a  signal  is  actually  present.  Again, 
using  the  asymptotic  approximation  and  the  same  nomenclature 


-CD  -o:> 


=  i  [1  -  $(d  -  a)j  [1  +  Ka)]  (29) 

Curves  of  versus  d  are  given  in  Fig.  3  for  various  values  of 
false-alarm  probability.  As  discussed  later  in  this  section,  these 
curves  must  be  straight  lines  on  normal  probability  scales.  The 
dotted  curves  are  shown  for  values  of  detection  probability  less  than 
50  percent,  to  indicate  that  the  asymptotic  decision  curve  approxi¬ 
mation  begins  to  br ^ak  down.  Actually,  an  error  check  was  made  for 
a  =  .01  and  P^  =  50  percent,  which  requires  d  =  2.6,  with  the  result 
that  the  probability  of  the  vector  (x^,  X2)  falling  in  the  error 
region.  Region  II,  is  1.506  percent.  Therefore,  the  correct  curve  of 
Pj  versus  d  would  curve  up  slightly  as  d  decreases  and  at  d  =  2.6 
would  show  P^  =  51.506  percent.  For  higher  values  of  detection 
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Fig. 3— Coherent  sinusoid;  signal  detectability  for  case  of  two 
possible  frequencies  (m  »  2) 


Detection  probability 
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probability  and/or  lower  values  of  false-alarm  probability,  this 
error  would  dei  .  very  rapidly. 

It  must  be  stressed  again  that  this  is  only  an  error  in  the 
analytical  results  to  facilitate,  a  closed- form  solution.  It  does  not 
reflect  an  error  made  by  the  system.  For  the  example  given 
above,  the  actual  system  has  a  detection  probability  1.5  percent  higher 
than  that  shown  in  Fig.  3.  A  negligible  error  i.3  made  in  the  calcu¬ 
lation  of  a,  since  the  density  function  as  seen  in  Fig.  2  is  situated 
so  far  from  Region  II . 

m-FREQUENCY  CASE 

If  the  number  of  possible  frequencies  at  which  the  signal  might 

appear  is  greater  than  2,  the  decision  surface  is  now  m-dimensional 

and  satisfies,  according  to  Eq.  (5),  for  equally  probable  frequencies 

m 

y  exp[log  L. ]  -  mk  (30) 

i=l 

where  k  is  the  preset  threshold,  and  log  is  the  quantity  described 
by  Eq.  (11)  for  each  frequency  In  order  to  calculate  the  per¬ 

formance  probabi' ities ,  it  is  necessary  to  find  the  joint  conditional 
probability  distributions  of  the  m  test  quantities  log  L^,,  i  =  1,  2, 

...,  m,  under  each  of  the  hypotheses  H  ,  H, ,  ,  ...,  H  ,  where  H, 

o  I  Z  m  1 

t  h 

indicates  the  assumption  that  a  signal  is  present  at  che  i  fre¬ 
quency  (1)^.  It  is  clear  from  the  analysis  of  Appendix  A  that  the 
result  will  be  an  m-diraensional  normal  distribution  of  uncorrelated 
variables,  and 

1.  if  no  signal  is  present,  all  of  the  m  quantities  log 

2 

i  =  1,  2,  ...,  m  will  have  mean  values  -d  /2; 
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2.  if  a  signal  is  present  at  uj.,  m  -  I  of  the  quantities  will 


2 

will  have  a  positive  mean  value  of  d  /2. 

The  m-dimensional  normal  distributions  must  be  integrated  over 

a  poition  of  the  decision  space  determined  by  the  uecision  surface 

of  Eq.  (30).  Although  the  actual  detector  will  form  the  sus'  of 

Eq.  (30)  and  compare  this  sum  with  the  threshold,  for  purposes  of 

analysis  it  is  necessary  to  make  the  approximation  to  the  asymptotic 

2 

decision  surface  as  was  done  in  the  m  =  2  case.  If  d  is  large  (due, 
for  instance,  to  a  long  observation  time),  then  the  mean  values  of 
the  quantities  mentioned  above  will  be  large,  either  negative  or 
positive.  In  Eq.  (30)  it  is  seen  that  if  one  term  is  larger  than 
the  others,  as  will  be  the  case  when  a  signal  is  present  at  one  of  the 
the  exponentiation  performed  will  markedly  accentuate  this  term 
and  als^o  depress  the  other  terms.  Thus,  the  detector  will  be  ef¬ 
fectively  making  its  decision  based  on  only  one  term.  This  is  equiva¬ 
lent  to  a  maximum  likelihood  detector  for  which  the  asymptotic  decision 
curves  are  appropriate. 

The  error  probabilities  may  easily  be  computed  by  combinatorial 

methods.  Since  the  m  quantities  log  i  =  1,  2,  ...,  m  are  sta- 

2 

tistically  independent,  and  for  large  d  the  asymptotic  decision 

curve  approximation  is  made  so  that  a  decision  results  according  to 

e.'’.ch  value  log  alone,  an  overall  false  alarm  will  result  if  any 

one  of  the  quantities  log  causes  a  false  alarm.  Each  test 

2 

quantity  is  normally  distributed  with  variance,  d  ,  and  with  mean 
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2  2 

values  -d  /2  or  +d'‘72  under  hypotheses  and  respectively.  Here 
it  is  understood  that  indicates  the  hypothesis  that  no  signal  is 
present  iti  the  i  quantity,  although  it  may  be  present  at  s’me  other 
frequency.  Thus,  for  each  test  quantity,  the  false-alarm  and  false- 
dismissal  probabilities  and  respectively)  are  given  by 


(31) 


where  the  asymptotic  threshold  value  log  mk  is  used.  The  detector 
rings  a  false  alarm  if  any  one  of  the  quantities  causes  a  false 
alarm.  The  probability  of  this  event  is  given  by 

a  =  1  ~  (1  -  0'^)"'  (32) 

for  small  <  <  1.  In  addition,  a  false  dismissal  occurs  wb^n 
all  of  the  m  quantities  do  not  exceed  the  threshold,  under  the  con¬ 
dition  that  there  is  a  signal  present  at  one  of  the  frequencies. 

Since  the  m  test  quantities  are  identically  distributed,  it  makes  no 
difference  at  which  one  the  signal  actually  appears,  that  is,  which 
of  the  quantities  log  produces  the  false  dismissal.  Hence,  a 
false  dismissal  occurs  whenever  there  are  m  -  1  correct  dismissals 
and  one  false  dismissal,  and  has  a  probability  of  occurring 

3  =  p.  (1  -  «  p. 

again  for  ma^  <  <  1. 


(33) 
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The  prcoabiiity  of  detec tiion,  ^  plotted  in  Fig.  4 

-2 

for  several  values  of  m  and  for  false-alarm  probabilities  of  10 

"  6 

and  10  .  It  is  seen  that  detectability  decreases  as  m  increases, 

which  is  reasonable  since  there  is  greater  chance  of  an  error  with 
a  large  number  of  independent  processors.  Also,  using  the  approxi¬ 
mations  of  Eqs.  (32)  and  (33),  Eq.  (31)  may  be  used  to  solve  for 
log  mk  in  terms  of  the  allowable  false-alarm  probability  Using 

this  value  of  threshold  to  compute  false-dism..ssal  probability  by 
Eq.  (31) 


1  -  3  =  -J-fi  +  $  d  -  f  1  -  —  ')  } 

21.  L  V  m  y  JJ 


(34) 


2  -1 

for  large  values  of  d  where  i  (•)  is  the  inverse  function  to  $(  ) . 

It  is  thus  seen  that  P,  *=  1  -  3  is  a  linear  function  of  d  on  normal 

d 

probability  scales,  since  the  argument  of  the  inverse  normal  proba¬ 


bility  function  is  constant  with  respect  to  d.  This  sheds  light  on 
the  linearity  of  the  detection  probability  curves  versus  d,  for 
constant  a-  Finally,  it  should  be  noted  that  the  loss  in  detecta¬ 
bility  resulting  from  doubling  the  number  of  possible  frequencies  m 
in  which  the  signal  may  appear  is  equivalent  to  subtracting  about 
.15  from  the  detection  index  at  very  low  values  of  ot,  and  a  slightly 
greater  amount  at  larger  values  of  a.  The  curve  for  m  =  1  is  the 
known  frequency  and  known  phase  detection  case. 


Miss  probability. 


Fig. 4— Coherent  sinusoid:  signal  detectability  for  case  of 
m  possible  frequencies 


Detection  probability 
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IV.  DETECTION  OF  AN  INCOHERENT  SINUSOID  OF  UNKNOWN  FREQUENCY 


The  more  realistic  case  of  the  detection  of  a  sine  wave  of  both 
unknown  frequency  and  phase  is  now  : -nsidered .  The  signal  again 
has  the  form  described  in  Eq.  (1),  hut  now  the  phase  (p  is  also  random 
and  uniformly  distributed  in  (0,  2n) .  The  resulting  likelihood 
ratio,  which  performs  the  average  over  both  frequency  and  phase, 
is  given  in  Eq.  (6).  The  quantity  in  the  summav.J  has  been  treated 
extensively  in  the  literature,''  ’  as  it  is  the  likelihood  ratio 
for  a  sine  wave  of  known  frequency  and  constant  amplitude  but  unknown 
phase.  For  white  gaussian  noise,  the  summand  becomes 


2tt 


p2TT 

dcp  L(v,uu^cp)  ^  C 


(35) 


o 

where  C  is  a  constant  proportional  to  the  output  signal-to-noise  power 

(4) 

ratio,  Iq(x)  is  the  zero-order  modified  Bessel  function,  '  and  the 
quantity  r^  is  given  by 


r, 

1 


2? 

N 

o 


\j ^  J  v(t)  cos  (Ju^t  dt  v(t)  sin  u)^t  dt  ^ 


Each  test  quantity  r^  is  formed  by  cross-correlating  the  re¬ 
ceived  signal  v(t)  with  the  two  quadrature -ccmponent  sinusoids  at 
the  i"  i-rcquency.  If  a  signal  at  is  indeed  present  in  v(t)  ,  it 

cannot  simultaneously  be  out  of  p*-  .se  with  both  of  the  quadrature- 
component  sinusoids;  consequently,  at  least  one  of  the  cross- 
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correlations  must  yield  a  non-zero  value-  r^  then  forms  the  normalized 
rms  value  of  these  two  cross-correlations-  The  quantity  r^^  may  also 
be  generated  by  passing  v(t)  through  a  filter  with  impulse  response 


h  ^(t) 
opt 


cos  U)^(T  -  t) 


0  <  t  <  T 
t  >  T 


(37) 


for  some  constant  K  followed  by  an  amplifier  with  a  gain  of  2P/N^K. 
This  may  be  seen  since  the  output  of  the  filter  at  time  T  is 
given  by  the  convolution  integral  as 

“.(T)  =  ¥-  f  v(T) 


N  J 
o  o 


cos  UJ.  T  dx 
1 


(38) 


The  final  likelihood  ratio  which  the  detector  must  form  becomes, 
for  equally  likely  frequencies 


m 


«V)  -  ^  c  I  I^(rp 


(39) 


1= . 


The  optimum  detector  consists  of  a  bank  of  narrow-band  f .iters,  each 
centered  at  one  of  the  m  possible  frequencies  which  the  signajL  might 
have.  Thus,  the  optimum  detector  forms  m  terms  r^,  i  =  1,  m 

according  to  Eq.  (36) ,  and  passes  each  r^  through  a  nonlinear  device 
that  takes  its  modified  Bessel  function.  These  quantities  are  then 
summed,  and  the  result  is  compared  to  a  threshold.  This  is  shown 
schematically  in  Fig.  5  on  the  following  page,  where,  for  simplicity. 


the  term  ~  has  been  removed  from  j£(y)  above  and  used  to  modify  a 


threshold  c.  The  final  threshold  is  called  k.  This  threshold  will  be 
adjusted  to  fix  the  false-alarm  probability  at  the  allowable  maximum, 
and  30  its  exact  form  here  is  not  important. 
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By  performing  the  analysis  first  for  the  simple  m  =  2  case,  the 


likelihood  ratio  of  Eq.  (39)  becomes,  with  the  constant  term 


Fig. 5— Form  of  optimum  detector 


The  detector  forms  this  quantity  and  compares  it  with  the  threshold 
k.  If  the  sum  exceeds  k,  the  signal-accept  decision  is  made.  If 
the  sum  does  not  exceed  k,  the  signal-dismiss  decision  is  made.  A 
decision  curve  may  be  formed  in  direct  analogy  to  the  coherent  de¬ 
tector  case,  and  for  m  =  2,  it  v;ill  satisfy 

Using  the  coordinates  r^^,  r^,  this  curve  is  plotted  in  Fig.  6 
for  various  values  of  k,  along  with  che  asymptotic  values 

r^,r2  =  (k  -  1)  (42) 

where  I  ^  (*)  is  the  inverse  modified  Bessel  function.  The  use  of 
o 

the  term  asymptote  is  retained  here  by  analogy  with  the  coherent  case, 
even  though  there  is  strictly  no  asymptotic  value  since  r^^  and  r^ 
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Fig. 6— Incoherent  sinusoid:  decision  curve  for  case  of 
two  possible  frequencies  (m  =  2) 
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are  always  non-negative.  By  the  asymptotic  values  here  Is  meant  the 
largest  value  of  r  under  the  decision  curve,  or  the  value  of  the 
decision  curve  when  one  envelope  is  zero  (note,  1^(0)  =  •  It  is 

seen  in  Fig.  6  that  the  decision  curve  does  not  appro  ch  its  asymptote 
as  quickly  as  in  the  last  case,  shown  in  Fig.  1.  However,  the  area 
of  Region  II  is  approximately  constant  for  different  values  of  k, 
and  thus  for  large  signal-to-noise  ratios  this  area  becomes  pro¬ 
portionately  less  important.  For  example,  later  calculations  will 

show  that  the  threshold  value  k  changes  very  rapidly  with  signal-to- 

2 

noise  ratio  and  false-alarm  probability.  As  d  (later  shown  to  be 
the  signal-to-noise  ratio)  increases,  or  a  decreases,  the  necessary 
value  of  k  increases  exponentially.  To  use  two  examples,  for  a 
detection  probability  of  70  percent  with  a  false-alarm  probability 
of  .01,  the  required  value  of  d  is  3.65,  and  k  =  17,000.  For  a  de¬ 
tection  probability  of  90  percent  with  a  false  alarm  of  .0001,  the 

•  8 

required  d  is  5.6,  and  k  =  5  x  10  .  For  such  large  values  of  k,  the 
error  in  using  the  asymptote  is  clearly  negligible. 

In  conjunction  with  the  examples  above,  it  should  be  remembered 
that  the  false-alarm  probability  a  occurs  in  each  observation  inter¬ 
val.  If  it  is  desired  to  calculate  from  this  probability  the  ex¬ 
pected  number  of  false  alarms  per  year,  a  is  simply  multiplied  by  the 
number  of  observation  intervals  in  one  year.  This  is  so  because  the 
occurrence  of  false  alarms  is  independent  from  one  interval  to  the 
next,  and  therefore  the  number  of  false  alarms  in  one  year  .^s  bi- 
noraially  distributed  with  mean  value  qo,  where  q  is  the  number  of 
observation  intervals  in  one  year.  For  example,  if  or  =  .0001  and 
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each  interval  lasts  5  min,  there  are  approximately  100,000  intervals 
per  year,  and  the  expected  number  of  false  alarms  per  year  is  10. 


PERFORMANCE  OF  THE  DETECTOR 


In  order  to  find  the  probabilities  of  false  alarm  and  detection, 

the  joint  distribution  of  the  test  quantities  r^,  i=l,  2,  ...,m 

mucL  again  be  found  under  the  various  hypotheses  H  ,  H, ,  . . . ,  H  as 

o  1  m 

defined  earlier.  Due  to  the  fact  that  the  correlators  of  Eq.  (36) 
process  the  received  waveform  for  only  a  time  T,  they  exhibit  an 
effective  bandwidth  of  approximately  2/T  cps.  More  precisely,  each 
correlator,  or  filter  and  amplifier,  given  by  Eq.  (37),  has  the  ef¬ 
fective  frequency  response  characteristic  [Ref.  5,  p.  314] 

‘.if) 


sin{  2rr  f 


jH.(f)  I  =  ^ - - 

"  ""o  2n|f-f.l? 


(f  >  0) 


(43) 


i'  2 


The  mirror  image  of  this  response  is,  of  course,  also  present  at 
negative  frequencies.  Thus,  if  two  possible  signal  frequencies  lie 
near  each  other,  there  will  be  some  overlap  in  their  correlator  re¬ 
sponses,  as  shown  in  Fig.  7.  It  is  seen,  however,  that  this  overlap 
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is  small  if  and  differ  by  more  than  2/T  cps.  When  T  is  large, 
tnerefore,  there  should  be  very  little  dependence  between  the  output 
signals  from  such  correlators.  This  is  shown  more  precisely  in 
Appendix  B,  where  the  joint  probability  density  function  of  the  outputs 
r^^  and  r^  is  calculated.  It  mav  be  seen  there,  for  example,  that  the 
correlation  between  r^  and  t2  is  given  very  closely  by 


<'u  "  “ 


2  -  mp  - 

<'“2  ■  “■!>  2 


(44) 


if  UJ2  is  close  to  ,  the  case  of  interest,  and  is  very  nearly  equal 
to 


T 

2  sin  uu  — 

'^^2  ^  - T —  ^2  ^  ^  ‘^l^ 

■^2  2 


if  and  uu^  are  widely  separated.  The  quantity  ^  measure  of 

the  dependence  given  by  Eq.  (B-21) ,  and  depends  on  the  cross  moments 
defined  in  Eqs.  (B-8)  and  (B-10' .  it  can  be  seen  that  two  conditions 
will  insure  that  this  dependence  is  negligible:  ja'^  -  j  T  >  >  1, 
and  for  each  frequency  o)^  T  >  >  1.  This  will  always  be  true  in 
practice.  Thus,  the  envelopes  are  independent,  and  their  joint 
distribution  under  each  hypothesis  can  be  obtained  from  Appendix  B. 
For  noise  only 


^1^2 


pCr^,r2/H  )  =  , 

d 


exp 


2d‘ 


(rj  +  r^)] 


(46) 


and  for  signal  plus  noise 


r  T 

12 


exp 


[■ 


2d 


/  2  ^  2  ^  ,4 

2  <*^1  ^"2 


>]  ^o( 


^2> 


pCrj^.r^/H^)  = 


(47) 
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2 

where  the  signal-to-noise  ratio,  or  detection  index  d  ,  is  identical 
to  that  of  the  coherent  case 


o 


Under  hypothesis  that  the  signal  appears  instead  at  oj^,  the  sub¬ 
scripts  in  Eq.  (47)  need  only  be  reversed. 

ACTUAL  FILTER  PERFORMANCE 

Having  obtained  a  set  of  correlators  as  the  desired  results  for 
the  optimum  detector,  it  is  instructive  at  this  point  to  digress  to 
the  subject  of  real  narrow-band  filters  in  order  to  compare  their 
performance  with  that  of  the  optimum  detector  above.  It  was  seen 
that  the  effective  bandwidth  of  the  correlator  acting  for  time  T  was 
approximately  2/T  cps.  This  is  the  minimum  attainable  bandwidth  for 
any  actual  filter  acting  only  for  T  sec,  as  is  well  known  from  the 
uncertainty  principle  [Ref.  1,  p.  21].  Thus,  if  an  actual  RLC  filter 
having  a  bandwidth  much  narrower  than  2/T  were  used,  the  filter  would 
be  expected  to  approximate  very  closely  the  optimum  detector.  The 
RLC  filter  is  considered  to  have  impulse  response  equal  to 

L  C 

h(T)  =  Kj^e  cos(a)^T  +  cp)  - — - !{ — t — ,  (49) 

v(t)  {  ?R  }  v^(t) 

2  2  2  2 

where  g  =  R/2L,  o)  =  o)  -  e  ,  yu  =  l/LC,  tan  p  =  e/uJ.  ,  and 

i  o  o  1 

=  2gu)^/^j^.  This  is  followed  by  an  amplifier  having  a  gain  of 
2P/N  ,  in  order  to  conform  with  Eq.  (38).  The  bandwidth  of  the 
filter  defined  by  the  half-power  frequencies  is  2e  cps.  At  time  T, 
the  output  voltage  is  given  by 
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2P 

V  (T)  =  “  v(t'‘  e 


-e(T-T) 


cos[a)^(T  -  t)  +  cp]  (It 


(50) 


It  is  shown  in  Appendix  C  that  for  an  input  noise  process  v(t)  with 

2 

spectral  density  of  level  N^/2  volts  /cps  over  a  band  (-C,  fO  that 
is  wide  with  respect  to  the  passband  of  the  filter,  the  variance  of 
the  output  process  (at  time  T)  is  given  by 


tfil  =■  ^  (eT  <  <  1)  (51) 

o 

The  autocorrelation  function  for  such  a  process  when  the  filter  has 
impulse  response  existing  only  for  time  T  is  also  given  there.  In 
order  to  find  the  signal-to-noise  ratio  at  the  output,  it  is  neces¬ 
sary  to  find  the  effect  of  the  filter  on  an  input  sine  wave.  As  shown 

2 

in  Appendix  B,  the  signal-to-noise  ratio  d  is  equal  to  the  square  of 

the  envelope  output  due  to  a  sine  wave  alone  at  time  T,  divided  by 

the  output  variance  due  to  the  noise  input.  (Due  to  the  non-gaussian 

statistics,  this  is  somewhat  different  from  the  coherent  case,  where 
2 

d  was  simply  the  square  of  the  shift  in  the  mean  value  of  the  output 
from  the  noise  only  to  the  signal  plus  noise  cases,  divided  by  the 
output  variance.)  The  envelope  output  may  be  found  as  shown  in 
Appendix  C.  This  value  is 


r(T) 


(eT  <  <  1) 


(52) 


o 

It  is  clear  that  when  this  is  squared  and  divided  by  the  value  in 

Eq.  (51),  the  same  value  is  obtained  for  the  signal-to-noise  ratio, 

2 

d  ,  as  before,  demonstrating  that  the  filter  performs  as  well  as  the 


optimum  detector  as  long  as  eT  <  <  i. 
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Consider  now  the  case  when  this  is  not  true.  In  practice,  with 
observation  times  in  the  order  of  hundreds  of  seconds,  it  is  very 
difficult  to  build  filter’s  with  bands  narrow  enough  to  satisfy  this 
condition.  Practicable  filters  will  have  too  wide  a  bandwidth,  which 
will  admit  too  much  noise  power,  thus  degrading  the  resulting  signal- 
to-noise  ratio.  This  is  shown  in  detail  in  Appendix  C,  where  using 
the  above  filter,  it  is  seen  that  the  variance  of  the  noise  output 
at  time  T  (provided  that  uj^T  >  >  1)  is 


t 


fil 


(53) 


when  eT  is  no  longer  much  less  than  unity.  This  is  seen  to  approach 
the  previous  value  as  eT  diminishes.  Also,  the  envelope,  output  at 
time  T  for  a  sine  wave  alone  is  shown  to  be 

r(I)  -  ^  [l  -  (54) 

o 

Squaring  this  and  dividing  as  before  by  the  variance  in  Eq.  (53) , 

2 

the  signal-to-noise  ratio,  or  detection  index,  for  the  filter 

detector  when  the  filter  has  a  bandwidth  no  longer  narrow  enough 
(eT  is  no  longer  necessarily  much  less  than  unity)  is 


P^T  2 


1  -  e 


■eT 


‘fil  N 


eT 


■eT 


(55) 


o  1  +  e 

It  is  thus  seen  that  a  degradation  factor  is  introduced,  the  ratio 
of  the  signal-to-noise  ratios  for  the  optimum  and  the  filter  de¬ 
tectors.  This  factor  has  the  value 

-eT 


d2  =  %  •  i 

cT 


1  +  e 


-eT 


(56) 
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and  it  seen  to  approach  unity  as  eT  approaches  zero,  and  it  approaches 
zero  as  eT  becomes  Large. 

2 

The  above  result  .  '‘r  may  be  used  directly  in  Eqs.  (46)  and 

(*+7)  ,  as  long  as  the  filter  outputs  are  approximately  independent. 

It  can  be  seen  from  the  derivation  in  Appendix  B  that  the  dependence 

2  2  2 

between  the  envelopes  is  given  by  the  quantity  where 

the  latter  quantities  are  cross  moments  between  the  quadrature  com¬ 
ponents  of  the  two  filters.  Using  the  technique  shown  there,  and  as¬ 
suming  that  two  identical  filters  are  used,  of  the  form  of  Eq.  (49),  it 
is  a  straightforward  although  tedious  matter  to  show  that  the  normal- 

ized  dependence  p'  =  given  by 

_  -4cT  _  -2eT  .  \  n,  ■  i 

_  ,2  e  -  2e  cos(u)„  -  ir  )  T  H*  1 

2  _  4e  _  „  ^  '  2  _ I'  _ 


4e^  +  (ui^  - 


1  -  e 


-2eT 


which  is  small  whenever  2e  <  <  1 5  that  is,  whenever  the  filter 

bandwidth  is  much  smaller  than  the  spacing  between  the  filter  center 
frequencies  (possible  signal  frequencies).  This  assures  approximate' 
independence  between  r^  and  and  Eq.  (55)  can  be  used  to  calculate 
the  signal-to-noise  ratio  when  real  RLC  filters  are  used  in  place  of 
the  optimum  correlators.  This  in  turn  may  be  used  to  calculate  the 
performance  probabilities  of  the  detector,  according  to  the  rel«^lons 


given  below. 


PERFORMANCE  PROBABILITIES 

Using  the  decision  curve  discussed  previously,  the  probabilities 
of  false  alarm  and  detection  may  now  be  calculated.  For  large  values 
of  signal-to-noise  ratio  at  the  output  of  the  detector  ,  the  as>’mptotic 
approximation  to  the  optimum  decision  curve  may  be  used,  as  the  type 
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of  argument  made  in  connection  with  the  coherent  detector  decision 

space  applicc  here  as  well.  That  is,  for  large  output  signal-to-noise 

ratios,  Pvegiou  II  of  Fig.  6  has  little  significance  since  the  centers 

of  the  distributions  in  question,  those  of  Eqs.  (46)  and  (47),  are 

located  far  from  this  legion.  In  particular,  if  there  is  no  signal 
til 

present  in  the  i  filter,  the  output  has  a  mean  value  r^  =  d.  If 

there  is  a  signal  present  in  this  filter,  the  output  has  a  mean  value 

2 

=  d  .  The  mean  values  are  separated  by  the  factor  d. 

In  the  actual  calculation,  the  false-alarm  probability  o'  is 
easily  computed  by  noting  the-  1  -  O'  is  the  probability  that  the  point 


lies  below  the  threshold,  given 


r  r  1  o  9  1 

O'  =  1  -  I  dr^  J  dr^  -J-  exp[-  (r^  +  r^J 
o  0  G  2d 


r  r 

iT' 

d^ 

2  c  .2 


(58) 


»  1  -  [l  -  exp  (-  -^  )  j 
L.  \  2d^  ' 


R# 


2d 


where 


b  =  l"^(k  -  1) 
o  " 


(59) 


(60) 


The  approximation  leading  to  Eq.  (59)  stems  from  the  fact  that 
o  is  required  to  be  small.  Using  Eq.  (59),  b  can  be  found  in  terms 
of  the  allowable  false-alarm  probability  per  observauion  time,  o*.  Thus 


Iog(2/Q*) 


(61) 


This  gives  the  threshold  value  k  of  Fig.  5 


k  =  1  4  I^j^V^d^  log(2/o*)  ] 


(62) 
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The  detection  piobability  can  similarly  be  found,  since  the 
miss  probability  (3=1-  P^)  is  Lae  probability  that  r^,  r^  lies 
below  the  threshold  if  a  signal  is  present,  say,  at  u)-. 


- 

•  2d^ 


*)] 


I  (r  ) 

o  2 


(63) 


exp 


where  Q(x,y)  is  the  Q-function,  given  by 


(64) 


Q(x,y)  =  J  t  e 
y 


x^) 


2 


I  (xt)  dt 
o  ' 


(65) 


which  has  been  tabulated  by  Marcum. 

A  useful  approximation  is  made  in  Eq.  (64)  by  neglecting 
the  exponential  term  as  being  much  smaller  than  one  and  then  using 
the  value  of  d  found  in  Eq.  (61).  This  gives  the  detection  proba¬ 
bility  in  terms  of  the  false-alarm  probability  and  the  detection 
index. 


Pd  =  Q  d,  si  1  log  2/a*  ) 


(66) 
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This  function  is  plotted  in  Fig.  8  for  several  values  of  allowable 
false-alarm  probability.  Also  presented  in  Fig.  8,  for  the  purpose 
of  comparison,  are  some  results  of  Fig.  3  for  the  coherent  case. 

For  large  values  of  signal-to-^noise  ratio,  v’here  t’r.e  curves  are  very 
accurate,'  the  incoherent  detector  yields  detection  probability 
curves  that  are  parallel  straight  lines.  This  indicates  that  the 
detection  probability  follows  a  gaussian  law,  which  may  be  seen  by 
using  the  common  approximation  to  the  Q- function  in  Eq.  (66),  given 
by  Rice  [Ref.  8,  p.  241].  It  is  important  to  realize  that  there  is 
a  constant  loss  between  the  coherent  and  incoherent  results.  In 
order  to  achieve  the  same  detection  probability  with  the  incoherent 
detector  as  with  the  coherent  detector,  the  detecti.cn  index  must  be 
increased  by  adding  about  d  =  .4.  For  example,  to  achieve  a  =  10  , 

=  .9,  the  coherent  detector  requires  d  =  5.15,  whereas  the  in¬ 
coherent  detector  requires  d  =  5.6.  Analysis  carried  out  in  terms 
of  known-phase  signals  would  finally  be  altered  by  adjusting  the 
required  d.  This  difference  could  presumably  be  made  up  by  in¬ 
creasing  observation  time.  The  attendant  small  increase  in  a  would 
have  little  effect.  This  amount  of  required  increase  is  found  to 
increase  slightly  with  higher  false-alarm  rates  (equivalent  to 
lower  thresholds) . 

This  result  indicates  that  in  the  case  of  detecting  a  sine  wave 
of  unknown  frequency  in  gaussian  noise,  the  problem  may  be  approached 
using  coherent  techniques,  which  are  more  convenient  to  work  with, 
and  then  introduce  some  constant  loss  factor  since  knowledge  of  the 
phase  is  not  available.  In  retrospect,  this  result  is  not  so 
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Detection  index,  d 

Fig. 8— Comparison  of  detectability  of  coherent  and  incoherent  sine  waves 
for  case  of  two  possible  frequencies  (m  =  2) 
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surprising  when  it  is  considered  how  the  detector  forms  the  optimum 
likelihood  ratio.  It  is  well  known  that  for  a  signal  of  known 
frequency  in  gaussian  noise,  the  necessary  increase  in  d  is  approxi¬ 
mately  .4  as  well  between  the  coherent  and  incoherent  detectors.  It 
seems  reasonable  that  the  same  sort  of  behavior  would  continue  to 
apply  when  a  set  of  such  variables  is  simply  added  in  each  case.  It 
is  useful,  however,  to  know  by  how  much  this  loss  factor  changes,  as 
given  above. 

DETECTOR  PERFORMANCE.  m-FREQUENCY  CASE 

If  the  signal  can  appear  at  any  one  of  m  possible  frequencies 
with,  say,  equal  probability,  the  detector  is  that  pictured  in  Fig.  5, 
and  the  decision  curve  satisfies 

m 

I 

i*l 

Again  assuming  that  |u)^  -  j  T  >  >  1  if  and  T  >  >  1  for 

all  i.j  so  that  the  m  quantities  I^(r^),  i  =  1,  2,  ...,  m  are  inde¬ 
pendent,  the  joint  conditional  density  function  of  the  m  quantities 
may  be  formed,  and  these  may  be  integrated  over  the  m-dimensional 
decision  space  to  obtain  the  false-alarm  and  detection  probabilities. 

In  a  manner  analogous  to  that  in  Section  III,  it  may  be  demonstrated 

2 

that  the  asymptotic  approximation  is  valid  for  large  values  of  d  , 
and  this  approximation  will  be  used  here.  Although  the  actual 
detector  does  form  the  likelihood  ratio  sum  and  compares  it 
with  the  threshold,  the  computation  uses  the  approximation  that 
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the  detector  forms  each  of  the  m  terms  r^  and  compares  each  of  them 
with  the  threshold  given  by  the  asymptote,  which  is 

r^  =  [k  -  (m  -  1)]  i  =  1,  2,  m  (68) 

Thus,  the  same  combinatorial  techniques  may  be  used  where  the  error 
probabilities  associated  with  each  quantity  r^  are  found.  These  then 
may  be  logically  combined  to  obtain  the  overall  error  probabilities. 
Each  quantity  r^  has  the  following  distribution 

P(r^/H^)  =  r./d^  exp(-r^/2d^)  (69) 

P(r./Hp  =  r./d^  exp[(-l/2d^)(r^  +  d^  ]  I^(r.)  (70) 

The  error  probabilities  may  be  formed  exactly  as  done  previously, 
with  the  result  for  each  test  quantity 


=  exp(-g^/2d^) 

(71) 

3.  =  1  -  Q(d,g/d) 

(72) 

where  g  =  ■  (m  -  1)].  Combining  these  for  the  m-frequency 

case,  the  overall  false-alarm  and  detection  probabilities  become 

a  =  1  -  (1  - 

=  m  exp[-g^/2d“]  (73) 

3  =  3^(1  -  «  3i 

so  that 

Pd  =  1  -  9i  =  Q(d.8/d)  (74) 
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Solving  Eqs.  (73)  and  (74)  for  in  terms  of  a*,  the  allowable  false- 
alarm  probability,  the  detection  probability 

Pj  =  Q  {i. -Jl  log  ^  )  (75) 

-2 

Curves  of  P,  versus  d  for  various  values  of  m  and  for  o'  =  10 
d 

“6 

and  10  are  given  in  Fig.  9,  along  with  some  of  the  coherent  de¬ 
tector  results  of  Fig.  4  for  comparison.  It  is  seen  thac  the  curves 
are  straight  lines,  indicating  that  the  detectability  follows  a 
gaussian  law  as  is  expected  from  the  normal  approximation  to  the 
Q-function.  Also,  it  is  seen  that  if  the  number  m  of  possible 
frequencies  at  which  the  signal  may  appear  is  doubled,  an  increase 
ill  d  of  between  .1  and  .2  is  required  in  order  to  achieve  the  same 
detection  probability.  This  very  slight  loss  is  encouraging  since 
in  practice  the  exact  frequency  would  not  be  known,  and  many  narrow- 
band  filters  would  have  to  be  built  in  order  to  cover  the  entire 
band  of  interest. 

The  results  seem  intuitively  satisfying,  since  in  Eq.  (75)  it 
is  seen  that  the  lack  of  knowledge  about  the  exact  signal  frequency 
requires  that  the  thresholds  be  set  slightly  higher  than  in  the  known 
frequency  case  (m  =  1).  Since  a  filter  is  used  for  each  possible 
frequency,  and  there  are  m  parallel  detectors  contributing  possible 
false  alarms,  a  higher  threshold  will  maintain  the  desired  false- 
alarm  rate.  Also,  the  appearance  of  the  simple  factor  d  in  Eq.  (75) 
indicates  that  if  a  signal  is  present  in  one  of  the  m  filters,  it 
alone  is  significant  in  causing  a  detection.  The  approximation 
of  Eq.  (74)  shows  this  by  neglecting  the  very  small  contribution  to 
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Fig. 9— Comparison  of  detectability  of  incoherent  sinusoid  for  case  of  m 
possible  frequencies  with  detectability  of  coherent  sinusoid  for 
case  of  two  possible  frequencies  (m  »  2) 
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the  optimum  detector,  by  passing  the 
function  device,  greatly  accentuates 
the  small  signals.  This  fact  causes 
present— to  predominate. 


filtered  signal  through  a  Bessel 
the  large  signals  and  depresses 
the  single  large  signal--if 
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V.  CONCLUSIONS 


The  results  of  the  above  analysis  indicate  that  it  is  not  neces¬ 
sary  to  work  in  terms  of  unknown  phase  detection  when  considering 
unknown  frequency  in  order  to  be  realistic.  For  the  case  of  a  finite 
number  of  known  frequencies  at  which  the  signal  may  appear,  it  is 
apparent  that  any  real  situation  may  be  approached  as  if  the  phase 
were  indeed  known  and  thus  work  with  gaussian  statistics.  This 
greatly  simplifies  the  analysis,  especially  for  a  large  number  of 
possible  frequencies.  An  adjustment  is  made  at  the  end  of  the  analysis 
to  account  for  the  actual  ignorance  of  the  signal  phase,  Ijt  this 
adjustment  is  quite  constant  for  many  different  situations,  thereby 
allowing  accurate  final  results. 

It  is  seen  that  the  optimum  detector  for  both  the  coherent  and 
incoherent  cases  is  a  nonlinear  device,  which  acutely  accentuates 
large  received  signals.  This  seems  quite  reasonable  under  the  cir¬ 
cumstances,  for  the  detector  is  testing  the  sum  of  a  large  number  of 
random  variables,  and  unless  the  detector  weighted  these  in  some 
manner  that  drew  attention  to  the  strongest  member,  or  equivalently 
suppressed  the  smaller  members,  the  false-alarm  rate  would  be  very 
high,  necessitating  a  higher  threshold  setting  which  would  reduce 
the  detection  capability. 

By  means  of  the  decision  plane  geometry  it  is  seen  that,  although 
the  optimum  detector  actually  tests  the  sum  of  the  filter  outputs 
against  a  threshold,  the  system  may  be  analyzed  as  if  each  filter 
output  were  separately  tested  against  the  threshold.  This  is,  of 
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course,  a  result  again  of  the  emphasis  placed  on  large  signals.  The 
system  then  becomes  one  of  a  band-splitting  type  which  divides  the 
frequency  band  of  concern  into  a  set  of  small  independent  bands  and 
interrogates  each  for  the  presence  of  a  signal.  The  approximation 
of  considering  the  optimum  detector  as  a  band-splitting  detector, 
which  is  made  clear  through  the  use  of  the  geometric  approach,  is 
very  g  for  large  output  signal-to-aoise  ratios.  It  is  seen  from 
the  d  jn-plane  representation  just  how  the  optimum  detector 

approaches  the  situation  of  m  approximately  orthogonal  signals  in 
noise.  The  asymptotic  approximation  finally  isolates  each  filter 
output,  permitting  m  independent  decisions. 

Even  though  the  power  levels  of  the  signal  and  noise  may  be  very 
unfavorable  at  the  hydrophone,  a  long  integration  time  will  increase 
the  signal-to-noise  ratio  to  a  useful  value.  This  requires  that  the 
noise  level  ai\d  signal  level  be  constant  for  the  duration  of  the 
observation  interval.  In  practice,  the  duration  of  each  observation 
interval  will  be  limited  by  the  noise  nor stationarity  and/or  the 
time  of  passage  of  the  target. 

(9) 

The  problem  of  unknown  noise  level  has  been  shown  by  Tuteur  ' 
to  be  of  critical  importance  in  an  actual  detection  situation.  He 
has  considered  the  effect  of  noise  uncertainty  on  the  detectability 
of  submarines  by  means  of  their  emission  of  broadband  signals  and 
has  demonstrated  that  unreasonable  detection  ranges  result  if  noise 
uncertainty  is  ignored.  He  suggests  a  long-term  average-power 
measuranent  of  the  noise  level,  which  would  include  and  presumably 
"swamp  out"  the  short-term  fluctuation  due  to  the  target's  own  wide- 


band  noise  as  it  passed  by  the  hydrophone,  as  a  means  of  estimating 
actual  noise  level.  This  technique  is  intuitively  appealing  in  the 
present  study  also,  since  the  sine  wave  power  will  have  little  effect 
on  the  long-term  wideband  noise  power  estimate.  There  is  need  fcr 
further  analysis  which  considers  the  detection  of  both  the  wide¬ 
band  noise  and  the  sine  wave  signals  simultaneously,  especially  in 
light  of  the  critical  nature  of  the  results  for  noise  level  uncertainty. 
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Appendix  A 

PROBABILITY  DISTRIBUTION  OF  THE  TEST  QUANTITIES 

It  is  required  to  find  the  joint  conditional  probability  density 
functions  for  the  test  quantities 


^1  "  \  " 

-  ^  s' 
2  -1 

K-l 

SjL  ^  V 

(A-l) 

^2  '  1-2  = 

1  , 

-1 

,  -1 

2  --2 

K 

(A- 2) 

where  the  quantities  are  as  defined  in  Section  III.  Since  the  only 
random  variable  appearing  in  these  expressions  is  the  vector  v,  which 
is  jointly  gaussian  (in  n-dimensions) ,  and  i ince  the  quantities 
Xj^,X2  are  formed  by  linear  superposition®  of  these  gaussian  variables, 
and  are  gaussian  random  variables,  rience  it  only  remains 
to  find  their  mean  and  variance  under  the  three  hypotheses,  H^, 
and  H^,  as  defined  in  Section  III. 

Let  the  notation 

(  )  (A-3) 

o 

indicate  the  statistical  average  of  the  quantity  in  brackets  under  the 
hypothesis  H^,  and  similarly  for  and  H2.  If  no  subscript  appears, 
then  V  has  already  been  changed  to  its  component  parts,  and  the  average 
is  that  over  the  uoise  portion.  It  is  clear  that  these  hypotheses 
imply 
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:  v(t)  =  n(t) 


:  v(t)  =  n(t)  -I-  s.j^(t) 


CA-4) 


H2  :  v(t)  =  n(t)  S^CO 

Thus,  using  standard  matrix  techniques,  the  conditional  averages 


\ 

='1  A 


K  ^  V 


-  2  +  ^1  ^  \  z  / 


(A-5) 


since  the  noise  has  zero  mean,  and  where,  as  shown  in  Section  III 


d  1  /  v-1  -  1  /  v’l  o 

—  =  2  -^1  -1  ■  2  -2  -  -2 


(A-6) 


d  ^  ,2  ^  ,  „-l  /  \ 

“  +  d+iLlK  (^n/ 


(A-7) 


\  "i 


)  -  -f  ) 


-  ^  +  s'  k'^  s 
2  -1  -  -2 


(A-8) 


d2 

-  T 
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In  Eq.  (A-8)  ,  ^  ”2  because  this  term  indicates 

the  cross-correlation  of  two  sine  waves  of  different  frequencies. 

In  particular,  using  the  conditions  of  Eq.  (16),  where  the  noise  is 
assumed  white,  there  results  (provided  that  |  ^1  "  ^^2^  T  >  >  1) 


n 


cos(UL)^t  +  cp^)  003(0)^^  + 


(A-9) 


^  0 

This  is  not  an  exact  relation,  except  when  |f^  -  f 2 |  is  a  multiple  of 
1/T  cps,  but  the  term  in  Eq.  (A-9)  does  not  increase  with  time  for 
reasonably  large  T,  whereas  the  term  -  -j”  grows  linearly  with  time. 
Thus,  the  cross-correlation  term  becomes  negligible  with  respect  to 
the  —  term  for  large  observation  times. 

To  compute  the  variances  of  x^  =  log 

var^  [X^l  =  (  (-  |-  +  k'^  v)  )  -  -i  d'* 

o'  '  o 
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varjj^rx^l  -  (-  f  +  s;  K-‘(s^  +  n)  ) 


■  4 


/  '  +  i-  . 
2  ^ 


i;  K-^  nf  )  - 


1  A^ 

4  ^ 


'  f  ^  i{  K‘^  ii  i;  K-‘  ( 


n  /  +  i;  (an'  )k'^  ii  -  f 


•'  v-1 

=  ll  K  s^ 


(A-11) 


=  d' 


Tlis  same  results  are  obtained  when  =  log  is  averaged  over  the 
hypotheses,  since  and  are  similarly  distributed.  It  still 
remains  to  find  the  correlation  between  Xj^  and  x^.  The  normalized 
correlation  coefficient  is 
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Rs  0 

Thus,  Xj^  and  X2  are  uncorrelated  and  consequently  independent.  The 
joint  distribution  can  then  be  expressed  according  to  the  results 
found  here,  and  these  are  given  in  Eqs.  (21),  (22)  and  (23). 
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Appendix  B 

JOINT  DENSITY  FUNCTION  OF  TH£  ENVELOPES  OF  TWO  CX)RI<i:LAIORS 


In  the  following  derivation,  the  input  is  assvimed  to  be  a 

2 

stationary  gaussian  noise  process  with  the  spectrum  level  N^/2  volts  / 
cps  and  band  limited  to  (-Q,  Cl) .  It  is  assumed  that  T  >  >  1/Q. 

Thus,  the  noise  may  be  considered  as  being  white,  and  the  following 
way  be  used  for  the  input  autocorrelation  function 

N 

\(t)  =  Y 


where  6(x)  is  defined  by 

OS 

J  f(x)  6(x  -  a)  dx  =  f(a)  (B-2) 

‘•06 


and  f(x)  is  any  function  continuous  at  x  =  a. 

As  discussed  in  Section  IV,  the  optimum  detector  must  form  the 
two  quantities  defined  by 


2p 

Ij^  =  —  J  v(t)  cos  u)j^t  dt 
”o  “o 

and 

^2  "  ^  f  V 

o  o 

for  the  optimum  detection  of  a  sine  wave  at  frequency  rad/sec. 
It  then  forms  the  envelope  ^  using  these  quantities  according  to 
Eq.  (36).  It  was  shown  there  that  the  correlations  performed  in 
Eq.  (B-3)  could  be  carried  out  by  a  "narrow-band"  filter  centered 


(B-3) 


at  u)^,  as  described  by  Eq.  (37).  In  the  following,  the  terms  "filter" 
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and  "correlator"  will  be  used  interchangeably,  referring  always  to  the 
device  that  performs  the  operations  of  Eq.  (36). 

Considering  further  a  second  filter  centered  at 


o  o 


(t)  cos  0)2*^ 


r 

^4  "  ?  Jo  V 


(B-4) 


Under  the  hypothesis  that  the  input  v(t)  consists  of  noise  alone, 
the  outputs  I^,  I2,  and  are  each  the  result  of  linear  orocessi. ' 
of  a  gaussian  process  and  hence  are  all  jointly  normal  r;*:  i<  v;  ables. 
It  is  only  necessary  to  fin.’  the  first  and  uecond  momintp  for  each 
in  order  to  specify  the  joint  distribution  completely. 

Each  quantity  has  zero  mean  since  n(t)  has  z<  -  mean.  The 
second  moments  are  calculated  as  follows 

~2  /  2P  - 

^I'viT'/J  ^  ri(a)  cos  cos 


^  t  da  ~  6(t  >  o)  cos 


coo 


T  COS 


■ff 

o  o 


dr  cos 


(B-5) 


P^T  ^ 
N 

o 


:>'-j 


The  result  is  not  exact,  but  as  T  increases,  the  rapidly  oscillating 
portion  of  the  integrand  contributes  very  little  to  the  variance. 
Similarly,  it  is  easily  shown  that 


-2  -2  -2  P^T  .  , 

^2  ^3  ^4 


The  covariances  are  calculated  in  a  similar  manner. 

^1^'’  ~  (  J  ‘^l^  '^l^ 


(B-6) 


^  1^’ 
“  N  J 


cos  sin  uj^T  dx 


(B-7) 


p2 

=  “  J  sin  2aj^T  dx  =  M- 


12 


o  o 

This  term  does  not  increase  with  time,  and  consequently  as  T  increases, 

it  v'ill  become  negligible  with  respect  to  the  variances  having  the 
2 

value  P  T/N  .  It  is  thus  set  to  zero  at  this  point. 


,2  ^T 


—  2p  r 

—  J  cos  OJ^X  cos  uu^x  dx 


o  o 


p2^  P  sin(a}2  *  ^  .  sinCo^  +  uj^^)  T 


N 


(uu^  -  T 


((ju^  +  T 


(B-8) 


p2^  sinCu)^  -  ujp  T  ^ 
n“  (uj  -  u>  )  T  "  ^13 

O  Z  i 
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The  approximation  in  Eq.  (B"8)  assumes  that  the  frequencies  vju^  and 
are  close  enough  together  so  chat  (uj^  +  j ,  thus 

making  the  sum  frequency  term  above  small  with  respect  to  the  dif¬ 
ference  frequency  term.  If  this  is  not  the  case  and  instead 

(^^2  ^  i  “2  ”  ‘“l  1  ’ 


p2^  2  sin  uj.jT 


(say  >  lUj^) 


(B-9) 


The  results  will  be  the  sane  in  either  case,  but  the  interesting 
case  cf  (1^2  close  to  will  be  considered  here,  and  Eq.  (B-8)  will 
be  used.  Similarly,  the  other  moments  are  easily  obtained. 


—  .  II  ^  p2  j*T 

^2^4  =  M.^3  ,  “  iT  J  ^  ° 

0  0 

(B-10) 

,  p2^  1  -  cosCuj^  -  T  _____ 

^l^v  “  T  (^2  -  (Jup  T  “  ^14  “  “^2^3 


These  quantities  define  a  symmetrical  covariance  matrix  M  like  that 
of  Rice  [Ref.  8,  p.  215],  which  may  be  inverted  to  form  the  infor¬ 
mation  matrix  M  ^ 


0 

"^13 

-^14 

0 

^14 

"^13 

■^13 

^^14 

'i' 

0 

■^14 

"^13 

0 

(B-11) 
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where 


t  = 


N 


and  A  =  \|t^  -  (m-^^ 


14' 


(B-12) 


=  '!» 


2 

'l2 


By  comparing  the  variances  and  the  cross  moments  defined  above,  it 
is  readily  seen  that  as  long  as  ju)^  “  j  T  >  >  1,  T  >  >  1,  and 
ua^  T  >  >  1,  all  of  the  cross  moments  will  become  negligible  with 
respect  to  \|f  for  large  T.  Thus,  M  ^  reduces  to  a  diagonal  matrix 

2 

indicating  that  the  random  variables  are  independent,  and  A  =  ^ 
results. 

Thus,  for  the  conditional  joint  density  function  under  the 
hypothesis  (under  which  all  variables  have  zero  mean) 


P(Il,l2’^3’^4'^o) 


/  2,2 
4tt  >j( 


iiT  +  ^^2  +  ^3  ^ 


(B-13) 


The  standard  transformation  may  be  made  to  the  Rayleigh  distribution 
for  the  envelopes  r^^  ^2  ^2  ^4  ?•  214], 

using 


1^  -  cos  cos  ©2 

I2  =  sin  -  r^  sin  G2 


(B-14) 


Using  standard  techniques  for  transformations  of  random  variables 
[Ref.  3,  p.  31],  the  following  is  obtained  for  the  conditional  Joint 
distribution  of  the  envelopes  r^^,  r^,  from  Eq.  (B-13) 


P(ri,r2/H^) 


(B-15) 
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cX'c  j  using  il(  of  Eq.  (B-12) 


2 

,2  P  T 
o 


(B-16) 


as  in  the  coherent  case.  Also,  this  form  for  r^^,  r^  is  seen  to  be 
exactly  that  for  the  r  of  Eq.  (36). 

Jk 

When  a  sine  wave  is  present  at  the  input  of  the  form  P  cosCou^t  +  Q 
for  any  angle  ",  it  is  clear  that  the  quantities  and  will  have 
their  means  shifted  from  zero  to  the  valuet 


13(1')  “  ^  cos(u)2t  +  §)  cos  oj^t 

o  “o 


dt  = 


p2^  P  sin(2u)2T  +  ?)  -  sin  §- 


cos  §  + 


(B-  17) 


d  cos  I  if  >  >  1 


X^CT)  =  ^  ^  cosCiu^t  +  5)  sin  u;  t  dt 

0  0 


p2j  r  cos  §  -  cos(2uj2T  + 

—  -  sin  5  +  — .  -  _ — - 

N  L  ^  2uj_T  J 

O  2 


(B-18) 


Ri  -d  sin  I 


if  :  >  1 


Due  to  signal  alone,  the  envelope  'will  have  a  value  of  apprcxi- 
2 

mately  d  .  If  the  above  mean  'elues  are  used  in  Eq.  (B-13)  and  the 
same  transfomation  as  in  E:/.  is  used,  the  joint  distribution 


of  r^,  may  be  easily  obtained  [Ref.  8,  pp.  214,  238].  Here  only 
the  final  result  is  given 


^l'^2 

p(rj^,r2/H2)  =  — 


-ij  (rj  +  +  db]  I  (r  )  (B-19) 

2d 


This  completes  the  derivation  of  the  envelope  statistics,  but  a  slight 

extension  is  now  made,  both  as  a  matter  of  interest  and  also  because 

it  yields  a  useful  measure  of  filter  output  dependence.  The  joint 

distribution  of  the  envelopes  r^  and  r2  when  the  dependence  is  not 

negligibly  small  is  given  here.  Thus,  starting  from  the  information 

matrix  of  Eq.  (B-11)  ,  retained  throughout,  letting 

only  and  be  zero,  which  is  reasonable  in  light  of  their 

definitions.  The  following  result  may  be  obtained  in  a  manner  similar 

to  Rice's  derivation  of  the  autocorrelation  function  for  a  filter 

output  [Ref.  8,  p.  216]  but  extends  the  analysis  to  the  case  when  a 

sine  wave  is  present  at  frequency  0^2*  dependence  between  r^^  and 

2  2 

r2  ia  seen  to  occur  as  ~  ►*]^3  ^]4’  ^2 

completely  independent  only  if  ^j^2  ” 

The  result  is  as  follows:  For  two  envelopes  r^^  and  r2  under 
hypothesis  H2.  assuming  only  that  y.j^2  ^34 

small 

p(r  r  /H  )  =  f  de  exp(  ^  )  exp[-  ^  P^e)] 

2nd  ff  o  2d  a 


(B-20) 
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where 


a  =  1 


(B-21) 


and 


4  2 


2r2  cos  6 


(B-22) 


A  closed-form  expression  for  this  was  not  found,  although  numerical 

techniques  could  be  used  to  evaluate  it. 

Under  the  hypothesis  of  noise  only,  the  only  change  is  that 
2  2  2 

3  (9)  =  r2/d  ,  which  permits  a  closed-form  evaluation 


P(ri 


“‘■’I 

d  CT 


[■ 


2d^a 


/  2  ^  2 
(ri  +  r^ 


>]'X 


Vl  -  o 


(B-23) 


which  agrees  v ’.th  Bendat  [Ref.  6,  p.  308]  and  provides  a  convenient 
check.  If  *12  =  '>•  then  a  =  1,  and  it  is  easily  seen  that  r^  and 
become  independent,  leading  to  Eqs.  (B-15)  and  (B-19)  under 
and  respectively. 
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Appendix  C 

AUTOCORREI.ATION  FUNCTION  FOR  THE  ACTUAL  NARROW- BAND  FILTER 


Consider  the  output  of  the  narrow-band  filter  with  impulse  re¬ 
sponse  as  in  Eq.  (49),  when  the  filter  is  connected  only  during  (0,T). 
The  scaling  factor  2P/N^Kj^  is  added  by  means  of  an  amplifier  in  order 
to  concur  with  Eq-  (50) 


h^(T) 


K.  e  ®\os(uu.  T  +  cp) 
0 


0  ^  T  ^  T 
Otherwise 


(C-1) 


The  input  process  to  the  filter  is  a  stationary  white  nc \se  process 

2 

of  the  level  N^/2  volts  /cps.  The  noise  could  as  well  be  band-limited 
to  a  band  (-Q,  Q) ;  for  as  long  as  the  passband  of  the  filter  is 
narrow  with  respect  to  this  band,  the  noise  could  still  be  assumed 
white. 

The  output  correlation  function  is  given  by  Lee  [Ref.  5,  p.  332] 


R(t)  =  J 


hj(u) 


du 


J 


R.(T 

1 


+  u 


a)  da 


(C-2) 


where  the  input  correlation  function  R.(t)  is  given  by  Eq.  (B-1)  and 
N  ^ 

equals  6(t).  Thus 


R(t)  =  (  )  J  e  ®^ccs(u)^v  +  cp) 

N  “o 
o 


e  cos(u)j^a  + 


9) 


N 

^  6(t  +  V  -  a)  da 


(C-3) 


_  ^  - 
“  N  J  ® 


e(2u-H-) 


cos(u}^v  +  9)  cos[u)^(t  +  u)  +  tp]  du 


(C-4) 


if  (t  +  u)  lies  in  (0,T) 
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and  equals  zero  if  t  -r  v  does  nit  lie  in  (0,  T) .  This  restriction 
on  T  +  u  requires  that  for  a  given  value  of  t,  the  range  of  v  is 
(”T,  T  -  t) ;  and  t  is  also  restricted  to  (0,  T) .  This  sets  the  upper 
limit  of  the  integral 


R(t)  “  ®  J  ®  '^^^[cos  tu.  T  +  cos(2to,  u  +  2cp  +  tu.  t)]  dv  (C-5) 

O  **0  111 

if  T  <  T,  and  zero  otherwise. 

The  rapidly  oscillating  portion  of  the  integrand  may  be  neglected 
as  it  contributes  very  little  to  R(t)  ,  provided  2ti^  }t  -  r  j  >  >  1.  Thus 


R(t) 


cos  UJ,  T 
1 


du 


(C-6) 


r 

N 


-eT 


cos 


U).  T 
1 


-2eT(T-T) 


] 


(C-7) 


and  the  variance  is  given  by  R(0)  as 

hn  •  ^  [1  -  (c-B, 

O 

For  a  sine  wave  input  to  the  circuit,  the  output  of  the  filter  at 
time  T  may  be  found  in  a  manner  similar  to  that  of  Eqs.  (B-17)  and 
(B-18)  of  Appendix  B.  The  result  is  that  the  envelope  has  the  value 


(for  2uj^T  >  >  I) 


(C-9) 


Consider  the  case  where  eT  <  <  1;  that  is,  where  the  filter  is 
very  narrow  band  relative  to  the  observation  time.  Then  Ec.  (C-8) 


reduces  to 
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(C-10) 


(C-11) 
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